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Abstract 

Current theories from biosocial (e.g.: the role of neurotransmitters in behavioral features), ecological (e.g.: cultural, political, 
and institutional conditions), and interpersonal (e.g.: attachment) perspectives have grounded interpersonal and romantic 
relationships in normative social experiences. However, these theories have not been developed to the point of providing a 
solid theoretical understanding of the dynamics present in interpersonal and romantic relationships, and integrative theories 
are still lacking. In this paper, mathematical models are use to investigate the dynamics of interpersonal and romantic 
relationships, which are examined via ordinary and stochastic differential equations, in order to provide insight into the 
behaviors of love. The analysis starts with a deterministic model and progresses to nonlinear stochastic models capturing 
the stochastic rates and factors (e.g.: ecological factors, such as historical, cultural and community conditions) that affect 
proximal experiences and shape the patterns of relationship. Numerical examples are given to illustrate various dynamics of 
interpersonal and romantic behaviors (with emphasis placed on sustained oscillations, and transitions between locally stable 
equilibria) that are observable in stochastic models (closely related to real interpersonal dynamics), but absent in deterministic 
models. 
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1 Introduction 

Interpersonal relationships appear in many contexts, such as in family, kinship, acquaintance, work, and clubs, to 
name a few. The manifestation of interpersonal relationships in society comes in many forms ranging from romantic, 
parent-child, friendships, comradeship, casual, friend-with-bencfits, soul-mates, dating to more recently Internet 
relationships. The most intriguing of all of these interpersonal relationships, which is also a dominant phenomenon 
and fundamental in human social life and interaction, is romantic relationship (PQ-0). 

Romantic relationships refer to the mutually ongoing interactions between two or more individuals. Recent works 
show that romantic relationships are more common among adolescents than has previously been assumed, with more 
than half of adolescents in the United States being involved in some form of romantic relationships ([3]- [2]). More 
than 70% of high school and college students report having had a special romantic relationship in the previous years, 
and also report more frequent interactions with romantic partners than with parents, siblings, and/or friends ( [4 - [5 j 
) . In the case of adult, the study of romantic behaviors may provide invaluable insight as to why majority of romantic 
relationships fail or do not make it to engagement and/or marriage ([3]). Surra and Hughes ([7]) found that more 
than half (54%) of couples in their studies exhibit unpredictable and nonlinear relational trajectories involving large 
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number of turning and tipping points. Partners identified events such as new rivals, unresolved differences, meeting 
partners family, and job changes as turning and tipping points that greatly changed and influenced the nature, 
quality and progress of the relationship. Similar findings by others scholars have appeared in various literatures 

m-m and m-m)- 

Research on romantic relationships among adolescents has gained traction with emphases on the quality of rela- 
tionships and their potential implications for positive and negative developmental and socio-psychological outcomes. 
Theories in biosocial (e.g.: effect of neurotransmitters in behavioral features) ([TO]- [H])i ecological (e.g.: cultural, 
political, and institutional conditions) ([32]), and interpersonal (e.g.: attachment) ([l3]-[14]) studies have grounded 
romantic relationships in normative social experiences and paradigms. However, these theories have not been de- 
veloped to the point of providing a solid theoretical understanding of the various dynamics present in romantic 
relationships; in addition, integrative theories are still lacking. The study of relationships has begun to hold both the 
artistic imaginations and interdisciplinary intellectual interests of various scholars in the fields of sociology, biology, 
neuroscience, psychology, anthropology, and mathematics ([TS]-[23]). Since experiments in these areas are difficult 
to design and may be constrained by ethical considerations, mathematical models can play a vital role in studying 
the dynamics of relationships and their behavioral features. However, there are few mathematical models capturing 
the various dynamics of romantic relationships. In this paper, we study both deterministic and stochastic models 
relationship from interpersonal perspective. 

Deterministic differential equations have been used extensively to study dynamic phenomena in a wide range of fields, 
ranging from physical, natural, biological to social sciences. The mathematical models capturing the dynamics of love 
between two people have recently gained attention among many researchers ( [20 - [28] ) who have provided extensions 
to Strogatz's seminal model. In an one-page influential work ([H]) and later in a book ([T3]), Strogatz applied a 
system of linear differential equations to study Shakespearean model of love affair of Romeo and Juliet. Rinaldi ([20 
and [21J ) , Sprott ( 23J), Liao and Ran ([24]) and Wauer et al ([28]) have investigated realistic perturbations and 
extensions of Strogatzian model by including features such as attraction factor ([20], [25], [53], and [55])) delay and 
nonlinear return functions ([55]), and three-body love affairs or love triangles ([53] and [55J). Rinaldi investigated 
the three mechanisms of love dynamics: instinct, return, and oblivion in ([22;), making the model more realistic 
due to the fact that it accounts for the growth of feeling from a state of indifference. In ([5U]), Rinaldi proposed 
a three dimensional model to describe the cyclical love dynamics of Laura and Patriarch and introduced nonlinear 
return and oblivion functions, and poems written by Patriarch are used to validate the dynamics. Gottman et al 
([57]) employed discrete dynamical models to describe the interaction between married couples; Liao and Ran ([541) 
studied time delays, nonlinear coupling and Hopf bifurcation conditions. Recently, Wauer et al ([55]) examined 
various models, starting with a time-invariant two dimensional linear and nonlinear models and concluding with 
time-dependent fluctuations in the source-terms and parameters. In previous papers ([5D]-[5T], and [53]-[55]), only 
dyadic interactions are considered, and other effects such as personalities and differential appeals of the individuals 
are ignored. As a result, learning and adaptation processes are ruled out. In this paper, we investigate stochastic 
dynamical models. But first, we summarized previous deterministic models ([15], |17j and |20j ) before developing an 
equivalent stochastic model. 

2 Models and Stability Analysis 

In this section, we study two models with two state variables. The variables Xi and X 2 are the measures of love 
of individual 1 and 2 for their respective partners, where positive and negative measures represent positive (e.g.: 
friendship, passionate, intimate) and negative (e.g.: antagonism and disdain) feelings, respectively. We first propose 
a deterministic system of differential equations to model the dynamic of romantic relationship, and later extend the 
model naturally to a stochastic dynamic model, where the deterministic rates become the stochastic rates. Using 
the typology of Strogatz ([T5]-[16]) and Sprott ([23]), the four romantic styles are summarized in figure 1. 

Summary of Figure 1: 

Region I: Eager Beaver: individual 1 is encouraged by his own feelings as well as that of individual 2 [oci > and 
A>0). 

Region II: Secure or Cautious lover: individual 1 retreats from his own feelings but is encouraged by that of indi- 
vidual 2 (at < and fa > 0). 

Region III: Hermit: individual 1 retreats from his own feelings and that of individual 2 (pti < and Pi < 0). 
Region IV: Narcissistic Nerd: individual 1 wants more of what he feels but retreat from the feelings of individual 



2 



II 



I 



a,. 



Ill 



IV 



Fig. 1. Typology and Characterization of Romantic Style 



2 (ai > and ft < 0). 

This classification allows us to characterize the dynamics exhibited by various combinations of a variety of romantic 
styles (see Fig. 1). Previous papers ([16]-[T7], [20] and [23]) have considered various dynamics using all possible 
combinations in the sign of parameters a, and ft. In this paper, we focus primarily on dynamics observed when 
individuals in region 77 interact with individuals in region III. We later focus on the features of stochastic model that 
are not present in the deterministic dynamics ([1Z])- For instance, Cherif ([E]) numerically showed that romantic 
relationships can exhibit exotic dynamics such as sustained oscillations while deterministic models show damped 
oscillations (e.g. stochastic resonance) and diffusion of trajectories between equilibria. In this section, we proceed 
as follows: deterministic (linear and nonlinear) models are provided, and are then followed by their stochastic 
equivalents. In the Deterministic Model section, we use standard deterministic linear romantic dynamic model ( |16j) 
and nonlinear model provided in ([17]). In Stochastic Dynamic Models, stochastic versions of the models studied in 
section 2.1 are investigated with special emphasis on behaviors not observed in the deterministic dynamics. We then 
provide an extension, where differential parameters are considered. 

2.1 Deterministic Models 

In ([17]), a more general model was proposed. In this paper, we focus on dyadic relational dynamics characterized 
by regions II and III interactions of figure 1 (secure or cautious love dynamics). To model the behavioral features 
of romantic dynamics, the following deterministic model is proposed: 



for (Xi,X 2 ) <E R x R, where cti > is non-negative, and ft and A^ i = 1,2 are real constant, respectively. These 
parameters are oblivion, reaction and attraction constants, respectively. For ft and Ai, we relax positivity condition. 
In the equations above (Eqs. 1- 2), we assume that feelings decay exponentially fast in the absence of partners. The 
parameters specify the romantic style of individuals 1 and 2. For instance, describes the extent to which individual 
i is encouraged by his/her own feeling. In other words, Q4 indicates the degree to which an individual has internalized 
a sense of his/her self-worth. In addition, it can be used as the level of anxiety and dependency on other's approval 
in romantic relationships. The parameters ft represent the extent to which individual i is encouraged by his/her 
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partner, and/or expects his/her partner to be supportive. It measures the tendency to seek or avoid closeness in a 
romantic relationship. Therefore, the term —aiXi say that the love measure of i, in the absence of the partner decay 
exponentially; and 1 / ai is the time required for love to decay. We propose this structure as opposed to model described 
in previous literature ([20 -[28 ), because romantic relationships (as in any interpersonal relationship) are not linear, 
especially return factors. The functional structure of return factors are motivated by dynamics often portrayed in 
romance novels and the tragic outcomes illustrated in them (Shakespeare's Romeo and Juliet, Patriarch's Canzoniere 
and Posteritati, Tolstoy's Anna Karenina, Flaubert's Madame Bovary, more recently the tragic myspace suicide of 
Megan Meier). The constant e in the return function can be interpreted as the compensatory constant. For example, 
the romantic dynamics between Laura Winslow and Steve Urkel in the sitcom Family Matter can serve as a popular 
media example. When Steve Urkel despairs, Laura Winslow feels sorry for him and her antagonism is overcome 
by feeling of pity. As a result, she reverses her reaction to passion. This behavioral characteristic is captured by 
the function of reaction or return function (e.g.: (3\X 2 (l — eX 2 ^j). This expression captures the compensation for 
antagonism with flattery, or pity, for positive and negative values of X 2 in fliX 2 (l — £X 2 ), respectively. For e = 0, 
the model reduces to the models proposed by Strogatz ([H]-[IS]), and others ([2DJ-IH] and [1S]-[2S]). In Strogatzian 
model of love affair, which corresponds to the case where e = in equations 1 and 2, the equilibrium point, (X\ , X 2 ) , 
is satisfied by the following equations ([29]): 



Xo = 



a 2 A 1 +/3 1 A 2 
a\a 2 - f3if3 2 
aiA 2 + fi 2 A x 
a\a 2 - f3iP 2 



(3) 
(4) 



For such a system (Strogatzian Model), we note the following theorem: 

Theorem 1 The equilibrium (Xi,X 2 ^J is non-negative and asymptotically stable if and only if: 

U\OL 2 



(5) 



and the equilibrium is otherwise unstable. The threshold Rd is the basic dyadic relationship threshold. It says that 
the system is stable if product of the ratios of reactiveness and oblivious coefficients is less than one. 

Proof. Since the system is linear and two dimensional, we can study its Jacobian. The Jacobian of the system (Eq. 
la-b) is given as: 



J = 



-a.\ Pi 
(3 2 -a 2 



(6) 



Since the trace r (J) = — (ai + a 2 ) < is non-positive, it suffices to show that the determinant A (J) is non-negative. 
The determinant A (J) condition is given as: 

A (J) = a x a 2 - j3x!3 2 > (7) 
From this, we obtain the necessary and sufficient condition for stability and positivity of the equilibrium: 



a\a 2 



(8) 



Equation 5 is equivalent to equation 8. Note that if Eq. 7 or equivalently equation 8 is not satisfied, the equilibrium 
point is a saddle point, which is unstable and the positivity of equilibrium point does not hold. □ 
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Fig. 2. Stability Domain. Stable dynamics are observed in regions A, B, and C, while unstable and unbounded trajectories 
reside in region D 

We can interpret theorem 1 as follows: for asymptotic stability, the squared geometric mean of the ratio of reactiveness 
to love and oblivion must be less than 1. Whenever this statement does not hold, Strogatzian model gives rise to 
unbounded feeling, which is obviously unrealistic. For that reason, we restrict our study to the stable condition and 
state the following corollary for the linear model. 

Corollary 2.1 For the Eq. 1, and e = 0, we have the following: 

(i) If Rd < 1 and P1P2 > or equivalently < Rd < 1, then the equilibrium point of the system does not admit stable 

focus point or center. The transients of Xj (t) cannot have damped oscillations or other cyclic dynamics, 
(ii) If Rd < I or Rd < and for some f3\p2 < 0, then the equilibrium point admits stable focus. 

The results of Collorary 2.1.1 can be illustrated by the figure 2, where a is the ratio of ct\ and ct2 (e.g. a — and R 
is the dyadic relationship threshold defined in equations 5. In region D, all equilibria are unstable and the dynamics 
give rise to unbounded (unrealistic) romantic feelings. Dynamics in the remaining regions (regions A-C) are stables. 
In region A and B, and C, we obtain stable focus and node, respectively. Later, we focus on the dynamics in region 
A since their stochastic counterparts exhibit different dynamics for some parameter ranges. Figure 3 shows the 
phase portrait and time series of romantic dynamics residing in region A in which the dynamics represent romantic 
behavior of secure and hermit individuals. 

For the nonlinear model where e 7^ 0, the collorary 2.1.1 does not hold entirely and theorem 1 has minor correction 
term. For more general models including the nonlinear case, the following statement holds. 

Theorem 2 The equilibrium (Xi,X2) is asymptotically stable if and only if: 



Rd 



P1P2 



did 2 < 1 



(9) 



aia 2 



where the correction term dj = 



with j — 1, 2 and g (it) is the linearized return function. The equilibrium is 
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Determinstic Model with p =5.7, a =+0.1, P 2 =-1, a 2 =+0.01 
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Fig. 3. shows a deterministic evolution of love measure for a romantic relationship between secured or cautious and hermit 
lovers. Top plot shows the phase portrait of equation 1-2 with e = 0, ax = 0.1, a.i = 0.01, f3i = 5.7, and 02 = —1, linearized 
near a stable equilibrium. The bottom panel shows the time series solutions of romantic feelings, which exhibit damped 
oscillations 

otherwise unstable. 

The proof of theorem 2 is similar to that of theorem 1. That is, we linearize around the steady state and change 
variables to obtain linear equations which then give rise to equation 9 (see [24 for a similar result). 

In the case of stochastic dynamics, our emphasis is placed on the dynamics that exhibit damped oscillations in 
the deterministic models (e.g. regions A and B), for the stochastic models fairly follow the deterministic behaviors 
of other equilibrium-type (e.g. stable nodes and limit cycle and/or centers) in region C. In the next section, we 
provide a stochastic dynamical system approach to illustrate some interesting dynamics observed in interpersonal 
and romantic relationships. 

2.2 Stochastic Dynamical Models 

In the previous mathematical papers on the subject, all factors of relationship are independent of each other and they 
consider time-invariant personalities and the appeal of individuals, ignoring long-term aging, learning, adaptation 
processes, fast fluctuation of feelings, and external forces and influences such as familial approval and disapproval of 
loved ones. The role of oxytocin or vasopressin in the behavioral features, cultural and institutional conditions, and 
attachment dynamics are also ignored. Accumulating all these forces as external factors that play major role on the 
quality of relationships, we consider stochastic variation of previous models ( 1 TJ ) to investigate the stochastic nature 
of romantic dynamics by considering the rates as stochastic rates. In the following section, we provide a method of 
using a deterministic formulation to derive a stochastic model. 

2.2.1 Derivation of Stochastic Love Dynamics 

In |17j . we outlined a similar method of deriving a stochastic equivalent of deterministic model. The method used 
to derive the stochastic differential equations for dynamical process naturally lead to Ito stochastic differential 
equations, as oppose to other stochastic calculi (e.g. Stratonovich) . This paper only summarizes the approximation 
procedures, and as before interested readers should consult the work of Kurtz ([3U]) for more detailed treatment 
of the diffusion equation approximation, which corresponds to a continuous time Markov process. We obtain the 
stochastic dynamical model for the processes by: 



(> 



Transition 


Rate 


X 1 -> Xj. - 1 


CK1.X1 


Xi -> JTi + 1 


/3iX 2 (1 - + 


X2 — ► X2 — 1 




X 2 -» x 2 + 1 


(1 - + A 2 



Table 1 
Transition rate 

(i) listing all the possible changes AX = [AXt, AX2] along with the probabilities for each change in a short time step 
At (see table 1); 



(ii) taking the expected changes E [AX] and covariance matrix E AX (AX) 



are calculated for the Markov process. 



Note that E [AX] (E [AX]) T = o (At 2 ) and can be ignored. The rates in table 1 become the conditional transition 

rates of the stochastic process, that is, P (Xtjt+At) — x i ~ = Xi) = — aiX\At+o (At) and so on. To each of the 
increments, we add and subtract its conditional expectation, conditioned on the value of the process at the beginning 
of the time increment of length At. This allows us to then decompose each increment into the sum of the expected 
value of the increment and sum of centered increment. That is, AXi = [— aiXi + P2X2 (l — £-^f) + ^1] At — AZ\ + 
AZ2 with the expected value of E (AXi) — \—oc\X\ + P2X2 (l — £X%) + Ai] At, where the centered increment 
AXi — E (AXi) is given as the difference of two increments, AZ 2 — AZ\. The terms AZi are the difference of two 
centered Poisson increments. These terms are then replaced by increment of Brownian motion dWi with corrected 
standard deviations or conditional variance. The stochastic equations of the process can then be expressed in a form 
easily comparable to their deterministic equation counterpart. 

Alternatively, we can also arrive at the stochastic model by dividing the expected changes and the square root of 
the covariance matrix by At. In the limit as At — > 0, the former becomes the drift term /i (t, Xi, X 2 ), and the latter 
becomes the diffusion coefficient D (i, X%, X2), respectively. Both procedures yield similar stochastic differential 
equations of the form: 



dX = n(t,X 1 ,X 2 )dt + D(t,Xx,X 2 )dW 



(10) 



where W = [Wi, ...,W^] T is an independent Wiener process. Notice that from the above formalism, the following 
statements are also true and can be verified: 



and 



E 



E 



OasAt -> 



E ( AX ^f )T ) - D (*, X U X 2 ) D (t, X X , X 2 f 



OasAt -> 



(11) 



(12) 



The alternative procedure relies on using the discrete deterministic model and using similar argument as in the first 
method. Using this approach, we arrive at a discrete stochastic model. Conditions (Eqs. 11-12) provide a justification 
for a weak approximation of moving from a discrete stochastic model to a continuous stochastic model. Eq. 



This weak approximation is equivalent to the convergence of a family of discrete state-space Markov chains to a 
continuous stochastic process. That is, for some class of smooth functions G : 5ft 2 — » 5ft and let the solution to the 
stochastic differential equations be X (T) at time T and the solution to the discrete stochastic equation be denoted 
by X A (T), then E [G (X (T))] — E[G (X A (T))] -> as At -> 0. This provides a definition for weak convergence 
of discrete to continuous stochastic differential equations. Kurtz ([30 ), and Kloden and Platen ([31]) have given 
detailed expositions on the methodology outlined in this section. 
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Fig. 4. shows the existence of sustained oscillation in the Stochastic model, while deterministic model does not exhibit such 
behaviors with the same parameter values e = 0, ai = 0.1, oli = 0.01, /3i = 5.7, and P2 = — 1. The deterministic system 
exhibits damped oscillations. It provides the contrast between the deterministic and stochastic models for the dynamics of 
romantic relationship. Left figure shows the distributions associated with the dynamics shown on the left. 

Using one of the methods sketched above, the stochastic differential equations describing the dynamics of romantic 
relationships are given as follows: 



dX 1 = [-a x X x + j3 x X 2 (1 - eXl) + A x ] dt - y /a x X x dW 1 + y / X X 2 (1 - eXj) + A x dW 2 (13) 
dX 2 = [-a 2 X 2 + fcXi (1 - eXf) + A 2 ] dt + JfeXi (1 - eXj) + A 2 dW 3 - ^Ja 2 X 2 dW 4 (14) 



for (X Xl X 2 ) £ R x R, where Wi,i — 1,...,4 are independent standard Wiener processes. This approach allows 
us to extend deterministic models to stochastic models. One can analyze the dynamics of the stochastic models 
with the help of the stability analysis of the deterministic equations. In fact, the solution to the deterministic 
model corresponds to the mean of the stochastic model. This framework provides a step towards understanding of 
the dynamics that are exhibited by the stochastic model. It should be noted that the dynamics of the stochastic 
differential equations 13-14 are closely related to the dynamics of the deterministic model (Eqs. 1- 2), but can exhibit 
important differences. Additional dynamics can emerge in the stochastic model for some parameter values (Figs. 
4-5) for some [3 X [3 2 < (e.g. for Strogatzian model). In (|17j). similar and more exotic behaviors were observed for 
nonlinear case and some of the results are included in herein, for both linear and nonlinear systems. 

For a stochastic Strogatzian love affair (for the case where e = 0) , the system exhibits sustained oscillations whereas a 
deterministic model shows damped oscillation (see figure 4). In the case of nonlinear model (e > 0) with appropriate 
conditions being satisfied, we observe various dynamics including those observed in the linear case (e.g. sustained 
oscillation via stochastic resonance). 

In figure 5, oscillations can persist or are sustained for some parameter values, whereas the deterministic equations 
exhibit damped oscillations (figures 4 and 5). Note that in figure 5, we observed both sustained oscillations at the 
equilibria and "jumping" or "switching" phenomena. In the case of "switching," there are trajectories which diffuse 
from equilibrium to equilibrium (see figure 5). A more detailed investigation of such dynamics (e.g. sustained and 
"jumping" oscillations or transition between locally stable equilibria), we hope, will be studied in the next paper, 
where a modification of multiple time-scale approach with Ito-Doeblin Formula (also known as Ito formula) can 
be employed. The system exhibiting sustained oscillations with transitions between local stable equilibria is only 
observed in fragile interpersonal and romantic relationships (figure 5), while robust interpersonal and romantic 
relationships shows only sustained oscillations (figure 4). These behavioral dynamics are dependent on the variance 
associated with the diffusion terms, which can easily be verified with multiple time-scale method outlined by Kuske et 
al ([SI]). We also observed (not show on the figures) that for some parameter values, the transitions (jump) between 
equilibria or sustained oscillations are transient and are not sustained for these values. In other cases, for different 
conditions, the trajectories visit most of the equilibria of the system in a sustained way. For example, when there are 
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Stochastic Model with B, =5.7, a,=+0.1, p 2 =-1 , a 2 =+0.1 
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Fig. 5. shows diffusion between two locally stable equilibria in stochastic dynamics of love affair, ft corresponds to the 
nonlinear stochastic dynamics with similar /3{ parameter values as in figure 3, where we use the following parameters: e 7^ 0, 
on = Q-2 — 0.1, Pi = 5.7, and @2 = — 1 • The distribution for X\ is bi-modal while that of X2 is unimodal. 

more than two locally stable equilibria satisfying corollary 2.1(H) condition, the trajectories can visit most (if not 
all) of the steady states in fragile relationships, hence exhibiting multi-modal distributions. Therefore, this paper 
illustrates the need for more mathematical analysis of interpersonal and romantic relationships from the perspective 
of nonlinear stochastic differential equations. 



2.2.2 Extension: Differential Romantic Style 



Using similar methods as outlined above, a more complex stochastic dynamic of love affairs can be proposed. The 



ansatz is as follows: 

dXi = [-aiXx + (3 X X<2 (1 - eX\) + alt ~ y/^X^dWi + JfkX 2 (1 - eX%) + A x dW 2 (15) 

dX 2 = [-a 2 X 2 + feXi (1 - eXf) + A 2 ] dt + [3 2 X X (1 - eXf) + A 2 dW 3 - ^a 2 X 2 dW 4 (16) 

da.\ — aio [an — a{\ dt + y/cti^dWs (17) 

da 2 = a 20 [a 2 \ - a 2 ] dt + ^/a 22 dW 6 (18) 

dpi = ho\fri-Pi]dt+y/?i^dW 7 (19) 

d(3 2 = f3 20 \pn - 2 ] dt + ^(h 2 dW & (20) 

(21) 



for (X\, X 2 , ai, a 2 , 0i, (3 2 ) £ R 2 x R 4 , where Wi,i — 1,...,8 are independent standard Wiener processes. Equations 
(15-21) represent a stochastic model of love affair with variability in romantic style which is impacted by environ- 
mental factors such as dislike of family members and/or friends against one's partner, ecological and institutional 
conditions, bio-sociological factors, etc... Variations due to these factors can also affect the measure of love in an 
unpredictable matter. One way to capture these variations is to have different variables for each one of the factors of 
interest. However, this will complicate the analysis, and render the analysis of the system intractable. The equations 
above simplify and assume that, for example, /3i(t,u>i, ...,u> n ), where u>i, i= l,..,n , are different romantic factors 
that affect the romantic style of individuals. Note that each of these parameter dynamics a±, a 2 , /3%, and (3 2 can be 
integrated and solved exactly. We can then take their limit cases (steady state conditions) and substitute them in 
equations 15-16 to reduce our system to equations 13-14. Assuming that romantic styles vary stochastically, one can 
also investigate the possible effect of the inclusion of these dynamics. One such effect is the spread in distribution of 
evolution of feelings. However, this analysis is left to the next paper and was not investigate herein. 
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3 Conclusion 



In this paper, we have considered both deterministic and stochastic models with nonlinear return functions. The 
stochastic model has a structure related to a deterministic model which allows us to study most of its dynamics 
through the lens of deterministic analysis. We have focused on particular subsets of interesting dynamics that are not 
observed in deterministic models. While a deterministic model exhibited damped oscillations with certain parameter 
values, the stochastic models showed sustained oscillations with the same parameter values. The results show that 
deterministic linear and nonlinear models tend to approach locally stable emotional behaviors. However, in the 
presence of stochasticity in the models, more complex and exotic patterns of emotional behaviors are observed. The 
stochastic differential equation extension provides insight into the dynamics of romantic relationships that are not 
captured by deterministic models, which assumes that love is scalar and individuals respond predictably to their 
feelings and that of others without external influences, such as ecological factors. Stochastic models capture the 
fluctuations due to the dynamics of love and that of external influences. This paper provides a new direction to the 
study of interpersonal relationship. The future direction toward more realistic mathematical and theoretical modeling 
of the dynamics of romantic relationship is possible through the lens of agent-based modeling, where community 
interaction is included, or integrative models (interpersonal, bio-sociological and ecological models are integrated) 
can be investigated. Calculation of the distribution of time indifference or apathy of individuals is a new possibility 
arising out of the stochastic dynamics of romantic relationships. The most fruitful direction from mathematical 
purview is developing methods to analysis systems that exhibit "stability boundary crossing" or "jump between 
locally stable equilibria" dynamics. The analysis of these dynamics, to the best of our knowledge, has not been 
investigated elsewhere and is worth studying in another paper. 
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